The significance of the imaginary part of the weak value 
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Unlike the real part ol the generalized weak value of an observable, which can in a restricted sense 
be operationally interpreted as an idealized conditioned average of that observable in the limit of 
zero measurement disturbance, the imaginary part of the generalized weak value does not provide 
information pertaining to the observable being measured. What it does provide is direct information 
about how the initial state would be unitarily disturbed by the observable operator. Specifically, we 
provide an operational interpretation for the imaginary part of the generalized weak value as the 
logarithmic directional derivative of the post-selection probability along the unitary flow generated 
by the action of the observable operator. To obtain this interpretation, we revisit the standard 
von Neumann measurement protocol for obtaining the real and imaginary parts of the weak value 
and solve it exactly for arbitrary initial states and post-selections using the quantum operations 
formalism, which allows us to understand in detail how each part of the generalized weak value 
arises in the linear response regime. We also provide exact treatments of qubit measurements and 
Gaussian detectors as illustrative special cases, and show that the measurement disturbance from a 
Gaussian detector is purely decohering in the Lindblad sense, which allows the shifts for a Gaussian 
detector to be completely understood for any coupling strength in terms of a single complex weak 
value that involves the decohered initial state. 

PACS numbers: 03.65.Ta,03.65.Ca,03.67.-a 



I. INTRODUCTION 

In their seminal Letter, Aharonov et al. [l[ claimed 
that they could consistently assign a particular value 
to an observable that was being weakly measured in a 
pre- and post-selected ensemble. To illustrate their tech- 
nique, they weakly coupled an observable A to a con- 
tinuous detector with an initial Gaussian wave-function. 
Normally, such a weak von Neumann coupling Q would 
approximately shift the mean of the Gaussian detector 
wave-function by the expectation value (V'll^lV'i) of ^ in 
the initial state IV'i), which would effectively measure A; 
however, they showed that by post-selecting a final state 
\'il)f) after the weak coupling, the mean of the Gaussian 
detector wave-function could be made to approximately 
shift by a complex quantity that they dubbed the weak 
value of the observable. 



(1) 



Notably, the weak value expression is not constrained 
to the eigenvalue range for the observable A, so it can 
become arbitrarily large for nearly orthogonal pre- and 
post-selections. 

This complex shift in the mean of the Gaussian de- 
tector wave-function was only approximate under weak 
von Neumann coupling and not directly observable, so 
its significance was not overtly clear; however, the Let- 
ter [l| also showed that both the real and imaginary 
parts of ([1]) could be operationally obtained from the 
linear response of the detector under separate conjugate 
observable measurements. The practical benefit of this 
observation was that one could amplify the response of 
the detector by making a clever choice of post-selection. 



which potentially allowed for the sensitive determination 
of other small parameters contributing to the evolution. 

After theoretical clarifications of the derivation in 0], 
experimental confirmation of such amplified detector re- 
sponse soon followed in optical systems [3, 0] ■ The am- 
plification has since been used successfully to sensitively 
measure a variety of phenomena to remarkable pre- 

cision, using both the real and imaginary parts of ([T]) as 
amplification parameters. Several theoretical extensions 
of the original derivation of the amplification [l^ - |27j and 
several proposals for other amplification measurements 
have also appeared [28[i32l |. In particular, it has been 
noted that how the amplification effect arises in such 
a continuous wave-function detector is not intrinsically 
quantum mechanical, but can also occur in classical wave 
mechanics fssj . which has prompted recent study into 
the mathematical phenomenon of superoscillations (e.g. 

MM)- 

Conceptually, however, the weak value expression ([1]) 
has remained quite controversial: since it is generally 
complex and not constrained to the spectrum of A, how 
should it be interpreted? Its primary interpretation in 
the literature has rested somewhat loosely upon the ob- 
servation that despite its anomalous behavior one can 
still decompose an expectation value through the in- 
sertion of the identity into an average of weak values, 

(ta.iai^,) = j:fm\^i)?m\A\A)/{^f\A)), which 

has the same form as decomposing a classical expectation 
value E(X\i) into an average of conditioned expectation 
values E{X\i) = Y.f PU\i)E[X\i, f). This observation, 
together with its approximate appearance operationally 
in weak conditioned measurements, make it tempting to 
interpret the weak value as a disturbance-free counter- 
factual conditioned average that can be assigned to the 
observable within the context of a pre- and post-selected 
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ensemble even when it is not strictly measured [36|-|39f . 

Supporting this point of view is the fact that when the 
real part of ([T|) is bounded by the eigenvalue range of A, 
it agrees with the classical conditioned expectation value 
for the observable [s^]- Moreover, even when the real 
part is outside the normal eigenvalue range, it still obeys 
a self-consistent logic [i^l and seems to indicate oddly 
sensible information regarding the operator A. As such, 
it has been used quite successfully to analyze and inter- 
pret many quantum-mechanical paradoxes both theoreti- 
cally and experimentally, such as tunneling time [4l| - |4^ , 
vacuum Cherenkov radiation [4^, cavity QED correla- 
tions Q , double-slit complementarity [47l . lisj , superlu- 
minal group velocities [i^, the N-box paradox [sO, [Ell 



phase singularities [52|, Hardy's paradox 
ton arrival time [57| . Bohmian tra ject ories 58l-l6l|. and 
Leggett-Garg inequality violations |62H64j . 

Arguably more important for its status as a quantity 
pertaining to the measurement of A, however, is the fact 
that the real part of ([1]) appears as a stable weak limit 
point for conditioned measurements even when the detec- 
tor is not a von Neumann-coupled continuous wave that 
can experience supcroscillatory interference (e.g. [63l - 
W\). As a result, we can infer that at least the real part 
of ([T]) must have some operational significance specifi- 
cally pertaining to the measurement of A that extends 
beyond the scope of the ori gina l derivation. This observa- 
tion prompted our Letter [68| showing that a principled 
treatment of a general conditioned average of an observ- 
able can in fact converge in the weak measurement limit 
to a generalized expression for the real part of ([T]) , 



Re^„ 



Tv{Pf{A,p.,}) 
2Tr(P/p,) 



(2) 



where {A^ pi\ = Api + piA is the anti- commutator be- 
tween the observable operator and an arbitrary initial 
state Pi represented by a density operator, and where 
Pf is an arbitrary post-selection represented by an ele- 
ment from a positive operator- valued measure (POVM). 
The general conditioned average converges to ^ pro- 
vided that the manner in which A is measured satis- 
fies reasonable sufficiency conditions [6^, [tOI that ensure 
that the disturbance intrinsic to the measurement pro- 
cess does not persist in the weak limit. 

It is in this precise restricted sense that we can oper- 
ationally interpret the real part of the weak value ^ as 
an idealized conditioned average of A in the limit of zero 
measurement disturbance. Since it is also the only ap- 
parent limiting value of the general conditioned average 
that no longer depends on how the measurement of A 
is being made, it is also distinguished as a measurement 
context-independent conditioned average. These observa- 
tions provide strong justification for the treatment of the 
real part of the weak value ([2]) as a form of value assign- 
ment [s^-IH, [tiI, [t^] for the observable A that depends 
only upon the preparation and post-selection [73j . 



However, we are still left with a mystery: what is the 
significance of the imaginary part of (HJ that appears 
in the von Neumann measurement, and how does it re- 
late to the operator Al We can find a partial answer to 
this question in existing literature (e.g. [l2], 113, EH, EH) 
that has associated the appearance of the imaginary part 
of ([ij in the response of the detector with the intrinsic 
disturbance, or back-action, of the measurement process. 
For example, regarding continuous von Neumann detec- 
tors Aharonov and Botero [13, P-8] note that "the imag- 
inary part of the complex weak value can be interpreted 
as a 'bias function' for the posterior sampling point [of 
the detector]." Furthermore, they note that "the weak 
value of an observable A is tied to the role of A as a 
generator for infinitesimal unitary transformations" (stI . 
p. 11]. Similarly, while discussing measurements of tun- 
neling time Steinberg [4l| states that the imaginary part 
is a "measure of the back-action on the particle due to 
the measurement interaction itself" and that the detector 
shift corresponding to the imaginary part "is sensitive to 
the details of the measurement apparatus (in particular, 
to the initial uncertainty in momentum) , unlike the [shift 
corresponding to the real part]." 

In this paper, we will augment these observations in 
the literature by providing a precise operational inter- 
pretation of the following generalized expression for the 
imaginary part of ([l}. 



Im^„ 



Tr(P;(-^[A,p,])) 
2Tv{Pfp,) ' 



(3) 



where [A, pi] = Api-- piA is the commutator between A 
and the initial state. We will sec that the imaginary part 
of the weak value does not pertain to the measurement of 
A as an observable. Instead, we will interpret it as half 
the logarithmic directional derivative of the post-selection 
probability along the flow generated by the unitary action 
of the operator A. As such, it provides an explicit mea- 
sure for the idealized disturbance that the coupling to A 
would have induced upon the initial state in the limit that 
the detector was not measured, which resembles the sug- 
gestion by Steinberg [4l|; however, we shall see that the 
measurement of the detector can strongly alter the state 
evolution away from that ideal. The explicit commutator 
in ([3]) also indicates that the imaginary part of the weak 
value involves the operator A in its role as a generator for 
unitary transformations as suggested by Aharonov and 
Botero [13, in contrast to the real part ^ that involves 
the operator A in its role as a measurable observable. 

To make it clear how the generalized weak value ex- 
pressions ([2]) and ([3|) and their interpretations arise 
within a traditional von Neumann detector, we will pro- 
vide an exact treatment of a von Neumann measurement 
using the formalism of quantum operations (e.g. f73 - [76j ) . 
In addition to augmenting existing derivations in the lit- 
erature that are concerned larg ely with understanding 
the detector response (e.g. [11 [ll [H, MM, M^), 
our exact approach serves to connect the standard treat- 
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ment of weak values to our more general contextual val- 
ues analysis that produces the real part [68l - [70| more 
explicitly. We also provide several examples that spe- 
cialize our exact solution to typically investigated cases: 
a particular momentum weak value, an arbitrary qubit 
observable measurement, and a Gaussian detector. As 
a consequence, we will show that the Gaussian detec- 
tor is notable since it induces measurement disturbance 
that purely decoheres the system state into the eigen- 
basis of A in the Lindblad sense with increasing mea- 
surement strength. Surprisingly, the pure decoherence 
allows the shifts in a Gaussian detector to be completely 
parametrized by a single complex weak value to all orders 
in the coupling strength, which allows those shifts to be 
completely understood using our interpretations of that 
weak value. 

The paper is organized as follows. In fjll]wc analyze the 
von Neumann measurement procedure in detail, starting 
with the traditional unconditioned analysis in ijll Al fol- 
lowed by an operational analysis of the unconditioned 
case in miB II and the conditioned case in ^11 B 21 After 
obtaining the exact solution for the von Neumann detec- 
tor response, we consider the weak measurement regime 
to linear order in the coupling strength in ijllll which clar- 
ifies the origins and interpretations of the expressions ^ 
and ([3]). We discuss the time-symmetric picture in ^IVI 
for completeness. After a brief Bohmian mechanics ex- 
ample in ijV Al that helps to illustrate our interpretation 
of the weak value, we provide the complete solutions for 
a qubit observable in ^VBI and a Gaussian detector in 
CI Finally, we present our conclusions in WII 

II. VON NEUMANN MEASUREMENT 

The traditional approach for obtaining a complex weak 
value [l| for a system observable is to post-select a weak 
Gaussian von Neumann measurement The real and 
imaginary parts of the complex weak value then appear 
as scaled shifts in the conditioned expectations of conju- 
gate detector observables to linear order in the coupling 
strength. To clarify how these shifts occur and how the 
weak value can be interpreted, we shall solve the von 
Neumann measurement model exactly in the presence of 
post-selection. 

A. Traditional Analysis 

A von Neumann measurement [11, Q unitarily couples 
an operator A on a system Hilbert space Jis to a momen- 
tum operator p on a continuous detector Hilbert space 
Hd via a time-dependent interaction Hamiltonian of the 
form, 

Hiit)^g{t)A®p. (4) 

The interaction profile g{t) is assumed to be a function 
that is only nonzero over some interaction time interval 



t £ [0, T]. The interaction is also assumed to be impulsive 
with respect to the natural evolution of the initial joint 
state p of the system and detector; i.e., the interaction 
Hamiltonian ^ acts as the total Hamiltonian during the 
entire interaction time interval. 
Solving the Schrodinger equation, 

ihdtU = HilJ, (5) 

with the initial condition Uq ~ 1 produces a unitary 
operator, 

Ut ^ cxp (^j^A (g) pj , (6) 

9= f dtg{t), (7) 
Jo 

that describes the full interaction over the time interval 
T. The constant g acts as an effective coupling parameter 
for the impulsive interaction. If the interaction is weakly 
coupled then g is sufficiently small so that Ut ~ 1 and 
the effect of the interaction will be approximately negli- 
gible; however, we will make no assumptions about the 
weakness of the coupling a priori. 

The unitary interaction ([6|) will entangle the system 
with the detector so that performing a direct measure- 
ment on the detector will lead to an indirect measure- 
ment being performed on the system. Specifically, we 
note that the position operator x of the detector satisfies 
the canonical commutation relation [x,p\ = ihl^, and 
thus will evolve in the Heisenberg picture of the interac- 
tion according to, 

{is®x)T ^uU^s®x)Ut, (8) 

= is®x + gA (g) id- 

As a result, measuring the mean of the detector position 
after the interaction {x)t = Tr((ls will produce, 

{x)t = {x)o+g{A)o. (9) 

Hence, the mean of the detector position will be shifted 
from its initial mean by the mean of the system observ- 
able A in the initial reduced system state, linearly scaled 
by the coupling strength g. For this reason wc say that 
directly measuring the average of the detector position x 
results in an indirect measurement of the average of the 
system observable A. 

The detector momentum p, on the other hand, does 
not evolve in the Heisenberg picture since [Ut, is'S^p] = 
0. Hence, we expect that measuring the average detector 
momentum will provide no information about the system 
observable A. 

As discussed in the introduction, however, when one 
conditions such a von Neumann measurement of the de- 
tector upon the outcome of a second measurement made 
only upon the system, then the conditioned average of 
both the position and the momentum of the detector can 
experience a shift. To see why this is so, we will find it 
useful to switch to the language of quantum operations 
(e.g. (t^-ItgI) in order to dissect the measurement in 
more detail. 
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Preparation 



Interaction 



Measurement 



the total probability density for detecting the position x, 



Ut 



t = 



• \x){x\ 

t = T 



FIG. 1. (color online) Schematic for a von Neumann mea- 
surement. An initially prepared system state pi and detector 
state |'!/'){^l become entangled with the von Neumann uni- 
tary interaction Ut (O over a time interval T. Measuring 
a particular detector position x after the interaction updates 
the detector state to |a;)(a::| and also updates the system state 
to Mx{pi), where Mx is an effective measurement op- 
eration that encodes the entanglement with and subsequent 
measurement of the detector. 



B. Quantum Operations 

1. Unconditioned Measurement 

As before, we will assume an impulsive interaction in 
what follows so that any natural time evolution in the 
joint system and detector state will be negligible on the 
time scale of the measurement. (For considerations of 
the detector dynamics, see [ill-) We will also assume for 
simplicity of discussion that the initial joint state of the 
system and detector before the interaction is a product 
state and that the detector state is pure, 



p = ® \ip){ipl 



(10) 



though we will be able to relax this assumption in our 
final results. Conceptually, this assumption states that a 
typical detector will be initially well-calibrated and un- 
correlated with the unknown system state that is being 
probed via the interaction. 

Evolving the initial state with the interaction unitary 
Ut © will entangle the system with the detector. Hence, 
subsequently measuring a particular detector position 
will be equivalent to performing an operation Aix upon 
the reduced system state, as iUustrated in Figure [U 



M,{pi) = Trd((i, ® \x){x\)UtpU^t) ^ M,P,Ml 



= {x\Uj 



(11) 
(12) 



where Trrf(-) is the partial trace over the detector Hilbert 
space, and is the Kraus operator associated with 
the operation Mx- Furthermore, since (xjV') = ^{x) is 
the initial detector position wave-function we find Mx = 
j da exp{—gadx)ip{x)\a){a\ ~ j daip{x — ga)\a){a\, or, 

more compactly, Mx = ipi^ — 5^)- 

If we do not perform a subsequent post-selection on 
the system state, then we trace out the system to find 



p{x) = Tr,(7W,(p,)) = Tr,(.E,p,), (13) 
{i,\ijl{\,®\x){x\)UT\i>). (14) 



Ex = MlMx 



where Trs(-) is the partial trace over the system Hilbert 
space. The probability operator Ex is a positive system 
operator that encodes the probability of measuring a par- 
ticular detector position x, and can also be written in 
terms of the initial detector position wave-function as 
Ex = \ip{x — gA)\'^. To conserve probability it satisfies 
the condition, J dx Ex = Is, making the operators Ex a 
positive operator-valued measure (POVM) on the system 
space. 

Consequently, averaging the position of the detector 
will effectively average a system observable with the ini- 
tial system state, 



/OO 
dxxp{x) =Trs{6p^), (15) 
'OO 

/OO 
dxxEx ^ {i'lUri'^s ® x)UT\i^), 
-OO 

= {x)ois +gA, 

where we see the Heisenberg evolved position operator 
dH]) naturally appear. 

Since the probability operators Ex are diagonal in the 
basis of A, then the effective system operator O will also 
be diagonal in the same basis. Hence, by modifying the 
values that we assign to the position measurements, we 
can arrange an indirect measurement of any system ob- 
servable spanned by {Ex} in the basis of A, including A 
itself, 



dx 



9 



(16) 



The chosen set of values (x— (a;)o)/.g are contextual values 
for A, which can be thought of as a generalized spectrum 
that relates A to the specific POVM {Ex} associated 
with the measurement context {A^^} |68l47C| . They are 
not the only values that we could assign to the position 
measurement in order to obtain the equality (jl6p . but 
they are arguably the simplest to obtain and compute, as 
well as the most frequently used in the literature. It is in 
this precise sense that we can say that the von Neumann 
coupling leads to an indirect measurement of the average 
of A in the absence of post-selection. 

The measurement of A comes at a cost, however, since 
the system state is necessarily disturbed by the opera- 
tions A^j; in order to obtain the probability operators Ex- 
The state may even be disturbed more than is strictly re- 
quired to make the measurement of A, which can be seen 
by rewriting the measurement operators in polar form, 
Mx — lJx\Ex\^^'^ , with the positive root of the probabil- 
ity operator |£J^|^/^ and an additional unitary operator 
Ux- This decomposition implies that A4x splits into an 
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effective composition of two distinct operations, 



(17a) 
(17b) 
(17c) 



We can interpret the operation £x that involves only the 
roots of the probability operator |£^a;|^/^ as the pure mea- 
surement operation producing E^. That is, it represents 
the minimum necessary disturbance that one must make 
to the initial state in order to extract a measurable prob- 
ability. The second operation Ux unitarily disturbs the 
initial state, but does not contribute to E^- Since only 
Ex can be used to infer information about A through the 
identity (jl6p , we conclude that the disturbance from lAx 
is superfluous. 

To identify the condition for eliminating lAx, we can 
rewrite the Kraus operator ([T^ using the polar form 
of the initial detector position wave- function tl){x) = 
eyiY>{i^s{x))^r{x), 



Mx 



exp(iV's(a; - gA))^pr{x - gA). 



(18) 



The phase factor becomes the unitary operator Ux = 
exp(i'0s(a; — gA)) for while the magnitude becomes 
the required positive root jiJ^:!^^^ = ijjr^x — gA) for Ex- 
Hence, to eliminate the superfluous operation Ux from a 
von Neumann measurement with coupling Hamiltonian 
([4]), one must use a purely real initial detector wave- 
function in position. 

For contrast, measuring only a particular detector mo- 
mentum p will be equivalent to performing a different 
operation Mp upon the reduced system state. 



A/-p(p,) = Trd((i, 



\p){p\)UTpUl) = Npp,Nl (19) 
'9P 



(pIV.). 



(20) 



The Kraus operator Np has a purely unitary factor con- 
taining A that will disturb the system, regardless of the 
form of the initial momentum wave- function (p| V') • More- 
over, the probability operator associated with the mo- 
mentum measurement has the form. 



± p — ±y j^i^p 



\{pW)?U 



(21) 



which can only be used to measure the identity Is. 

For completeness we also briefly note that the conju- 
gate Kraus operators Mx and Np are related through a 
Fourier transform, 

1 

Np^—= dxe-'P'-'/^Mx, (22a) 
V27rft j-oo 

1 f°° 

Mx^-= dpe^P^'/'^Np, (22b) 



• • 

Pi • *^Vf{Mxi^)) 

Mx 



t = T T2 



FIG. 2. (color online) Schematic for a sequence of two indi- 
rect measurements. After the von Neumann interaction and 
measurement of x illustrated in Figure [T] that produces the 
effective measurement operation A4x upon the initial system 
state, a second detector interacts impulsively with the sys- 
tem for a time interval T2 — T. The second detector is then 
measured to have a particular outcome /, which updates the 
system state to 'Pf{A4x{pi)), where Vf is another measure- 
ment operation. Taking the trace of the final system state 
will then produce the joint probability densities ((24)). 



and that both detector probability operators can be ob- 
tained as marginals of a Wigner quasi-probability opera- 
tor on the system Hilbert space, 



1 f°° 

Wx,p^— dye^^Py^^Ml^^M, 
dpWx,p, 



Ex 



dx ^Vx^p' 



x—y, (23a) 
(23b) 
(23c) 



In the absence of interaction, then the Wigner quasi- 
probability operator reduces to the Wigner quasi- 
probability distribution W(x,p) for the initial detector 

state, Wx,p W{x,p)is. 



2. Conditioned Measurement 

To post-select the system, an experimenter must per- 
form a second measurement after the von Neumann mea- 
surement and filter the two-measurement event space 
based on the outcomes for the second measurement. In 
other words, the experimenter keeps only those pairs of 
outcomes for which the second outcome satisfies some 
constraint. The remaining measurement pairs can then 
be averaged to produce conditioned averages of the first 
measurement. 

If we represent the second measurement as a set of 
probability operators {Pf} indexed by some parameter 
/ that can be derived analogously to from a set of 
operations {V/} as illustrated in Figure [5J then the total 
joint probability densities for the ordered sequences of 
measurement outcomes (x, /) and (p, /) will be, 

pixj) ^ TTsiPfMxip,)) = TisiExjM^ (24a) 
Pip, f) = Trs(P/A/'p(p,)) = Trs(FpjpO, (24b) 
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where the joint probability operators, 



where, 



Fpj = ATtp^TVp. 



(25a) 
(25b) 



are not simple products of the post-selection Pf and the 
probability operators or ([?!]) . Those operators can 
be recovered, however, by marginalizing over the index /, 
since the post-selection probability operators must satisfy 
a POVM condition J2fPf = ^s- 

The joint probabilities (j24p will contain information 
not only about the first measurement and the initial sys- 
tem state, but also about the second measurement and 
any disturbance to the initial state that occurred due to 
the first measurement. In particular, the joint probabil- 
ity operators can no longer satisfy the identity 
due to the second measurement, so averaging the proba- 
bilities (|24p must reveal more information about the mea- 
surement process than can be obtained solely from the 
operator A, the initial state pi, and the post-selection 
Pf. As a poignant example, the unitary disturbance Ux 
in (|17p that did not contribute to the operator identity 
(1161) will contribute to the joint probability operators, 
E,j = \E,\^/^U^,PfU,\Ex\^/\ 

The total probability for obtaining the post-selection 
outcome / can be obtained by marginalizing over either 
X or p in the joint probabilities. 



p{f)= I dxpixj)= I dppipj), (26) 

(27) 



= TTs{Pf£{pi)), 
£{p,)^^Td{^JT{p^<E)\i^){^\)Ul), 



where the operation £ is the total non-selective measure- 
ment that has been performed on pi. Since £ is not the 
identity operation, it represents the total disturbance in- 
trinsic to the measurement process. It includes unitary 
evolution of the reduced system state due to the inter- 
action Hamiltonian Q , as well as decoherence stemming 
from entanglement with the measured detector. 

By experimentally filtering the event pairs to keep only 
a particular outcome / of the second measurement, an 
experimenter can obtain the conditional probabilities. 



p{p\f) 



pjxj) _ Tr,(P/Al.(p,)) 

P{f) Ti- s{Pf£{M) 

p{pj) _ TrsiPfMpiM) 



Pif) 



Tr,(P/f(pO) 



(28a) 
(28b) 



which can then be averaged to find the exact conditioned 
averages for the detector position and momentum. 



f{x)T 
f{p)T 



d,.pi,\f)^TlAPfMEl}l^ (29a) 
Tr,(P,£(p,)) 

rfl^MpI/)-^^^^^^, (29b) 
^^^Vi^l^^ Tr,(P/£(p,)) 



XriM = Tr<j((i, (g, x)UT{p^ |^)(V'|)C/^), 



(30a) 
(30b) 



are detector averaging operations that affect the system 
state before the measurement of the post-selection is per- 
formed. It is worth noting at this point that we can relax 
the assumption (|10p made about the initial state in the 
exact operational expressions (P7|) and (PH)) . Similarly, if 
different contextual values are used to average the condi- 
tional probabilities in ([29]) . then corresponding detector 
observables with the same spectra will appear in the oper- 
ations pop in place of x or p; for example, averaging the 
values a{x) = {x— {x)o)/g used in (|16p will replace the de- 
tector observable x in ([5D|) with a. = dxa{x) |a;)(a;|. 

To better interpret pO|) . we bring the detector opera- 
tors inside the unitary operators in pO|) using the canon- 
ical commutation relations as in ([8l). 



^t(A) 



Xip^) - 



g£{{A,p,}l2), 



(31a) 
(31b) 



which splits the Xt operation into two operations but 
only changes the form of Vt- The operation proportional 
to g disturbs the symmetrized product {A,pi}/2 = 
{Api -\- Pi A) 12 of the initial system state with the oper- 
ator A, while the operations, 

X{pi) = TYd{lJT{p^ ® {x, |V^)(^|}/2)C/;), (32a) 
V{p,) = Tra{UT{p^ ® {p, K') (V'|}/2)f/;), (32b) 

disturb the symmetrized products of the initial detector 
state with the detector operators. 

The form of the equations (|3ip clearly illustrates how 
the post-selection will affect the measurement. If the 
post-selection is the identity operator, P/ = Is, then the 
unitary operators Ut causing the total disturbance of 
the initial state will cancel through the cyclic property of 
the total trace in (|29p . leaving the averages in the initial 
states that were previously obtained. 



{x)t = (x)o + .g(A)o, 
ip)T = {p)o- 



(33a) 
(33b) 



In this sense, commuting the detector operators x and 
p in ((30)) through the unitary operators to arrive at (|3T|) 
is equivalent to evolving them in the Heisenberg picture 
back from the time of measurement T to the initial time 
in order to compare them with the initial states. How- 
ever, the presence of the post-selection operator Pf will 
now generally spoil the cancelation of the unitary oper- 
ators that is implicit in the Heisenberg picture, leading 
to corrections from the disturbance between the pre- and 
post-selection. 

The symmetrized products in pip indicate the mea- 
surement being made on the initial states of the sys- 
tem and detector, which is then further disturbed by 
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the unitary operators TJt as a consequence of the cou- 
pling Hamiltonian The post-selection both condi- 
tions those measurements and reveals the disturbance, 
which corrects each term in (j33|) . yielding the final exact 
expressions, 



Tr,(P;A-(p,)) 



(34a) 



2Tr,(Py£(pO) 



(34b) 



III. THE WEAK VALUE 

If it were possible to leave the system state undisturbed 
while still allowing the measurement of A, then we would 
naively expect the disturbance £ to reduce to the iden- 
tity operation. Similarly, we would naively expect the 
operations X and V would reduce to {x)o and (p)o mul- 
tiplying the identity operation, respectively. As a result, 
the conditioned averages p4p would differ from the un- 
conditioned averages (p3|) solely by the replacement of 
the average {A)q with the real part ([2|) of the complex 
generalized weak value expression. 



-Am — 



TTsjPfAp,) 

'^^s{PfP^) ' 



(35) 



Since this expression depends solely upon the initial state 
Pi, the post-selection P/, and the operator A, we are 
naively tempted to give ReAu, an intuitive interpretation 
as the ideal conditioned expectation of A in a pre- and 
post-selected state with no intermediate measurement 
disturbance. However, it is strictly impossible to remove 
the disturbance from the measurement while still mak- 
ing the measurement, so we cannot rely on this sort of 
reasoning. We can make a similar interpretation in a re- 
stricted sense, however, by making the coupling strength 
g sufficiently small to reduce the disturbance to a minimal 
amount that still allows the measurement to be made. 

To see how the operations £, X, and V in (|27l) and 
P2p depend on the coupling strength we expand them 
pcrturbatively. 



n! \ih 



(36a) 



^(pO = E3f#)"^^A^(adA)"(pO, (36b) 



! \ih/ 2 

n—u 

^('^^) = E (1)" (p"+')o(adi)"(A), (36c) 

n=0 

where the operation (adA)(-) = [A, ■] is the left action of 
A in the adjoint representation of its Lie algebra, which 
takes the form of a commutator. That is, (adA) explic- 
itly describes how A disturbs the initial state due to the 
interaction that measures it. 



The initial detector state plays a critical role in ((36)) 
by determining the various moments, (p")o, (p"'*^^)o and 
{{p^,x}/2)o that appear in the series expansions. No- 
tably, if we make the initial detector wave- function purely 
real so that it minimally disturbs the system state then 
all moments containing odd powers of p will vanish. We 
conclude that those moments of the disturbance opera- 
tions are supcrffuous for obtaining the measurable prob- 
abilities that allow the measurement of A, while the mo- 
ments with even powers of p are necessary. 

After expanding the corrections ([M)) to first order in g, 
we obtain the linear response of the conditioned detector 
means due to the interaction, 

^(.)o + A (fa;})oT^-(f/(-^dA)(^-)) (37a) 



9- 



ih 2 Tr,(P/p, 
Tr,{Pf{A,p,}) 
2Tr,(P/p,) ' 



Measurements for which this linear response is a good 
approximation are known as weak measurements. 

After introducing the complex generalized weak value 
psp , we can write the linear response formulas in a more 
compact form, 

f{x)T = (x)o + l^^^y^ i2lmA^)+gReA^, (38a) 



/(p)t = (p)o + |(p'>o(2W^), 



(38b) 



in terms of not only its real part, but also twice its imag- 
inary part. 

If the initial detector position wave-function ip{^) is 
purely real, so that the measurement is minimally dis- 
turbing, then {{p, x}/2)q will vanish, leaving only ReA^ 
in f{x)T as we naively reasoned before. However, the 
term proportional to 2Im^„, will not vanish in f{p)T to 
linear order in g, making it an element of measurement 
disturbance that persists even for minimally disturbing 
weak measurements. 

These linear response formulas for the von Neumann 
measurement have also been obtained and discussed in 
the literature with varying de gree s of generality and rigor 
(e.g. 0,i[Ii,[il[Ii,|iiMliI|23i43)- However, our 
derivation has a conceptual advantage in that we see ex- 
plicitly how the origins of the real and imaginary parts 
of the weak value differ with respect to the measurement 
of A. We are therefore in a position to give concrete 
interpretations for each part. 

The real part ^ of the weak value Re^^ stems di- 
rectly from the part of the conditioned shift of the de- 
tector pointer that corresponds to the measurement of A 
and does not contain any further perturbation induced 
by the measurement coupling that would be indicated 
by factors of (ad A). As a result, it can be interpreted 
as an idealized limit point for the average of A in the 
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initial state pi that has been conditioned on the post- 
selection Pf without any appreciable intermediate mea- 
surement disturbance. To support this point of view, we 
have also shown in j68l - l70j that ReA^ appears naturally 
as such a limit point for minimally disturbing measure- 
ments that are not of von Neumann type, provided that 
those measurements satisfy reasonable sufficiency condi- 
tions regarding the measurability of A. 

The imaginary part ([3]) of the weak value ImA^, on 
the other hand, stems directly from the disturbance of 
the measurement and explicitly contains (adA), which 
is the action of A as a generator for unitary evolution 
due to the specific Hamiltonian (U). The factor 2ImA^ 
appears in (j38p along with information about the initial 
detector momentum that is being coupled to A in the 
Hamiltonian (|4]) , as well as factors of fi, in stark contrast 
to the real part. How then can it be interpreted? 

The significance of 21uiAw becomes more clear once we 
identity the directional derivative operation that appears 
in its numerator, 

5a{-) = -i(adA)(.). (39) 

That is, (5yi(pi) indicates the rate of change of the initial 
state pi along a flow in state-space generated by A. 

The directional derivative should be familiar from 
the Schrodinger equation written in the form dtp — 
\H,p\/ih = Sq{p), where the scaled Hamiltonian Ct = 
H /his a characteristic frequency operator. The integra- 
tion of this equation is a unitary operation in exponential 
form p{t) = exp(tfc)(/5(0)) = exp(— i<j2)p(0) exp(iir2) 
that specifies a flow in state space, which is a collection 
of curves that is parametrized both by a time parameter 
t and by the initial condition p(0). Specifying the initial 
condition /d(0) = p^, picks out the specific curve from the 
fiow that contains pi. The directional derivative of the 
initial state along that specific curve is then defined in 
the standard way, dtp{t)\t=o = Sn{pi). 

The fact that the quantum state space is always a con- 
tinuous manifold of states allows such a fiow to be de- 
fined in a similar fashion using any Hermitian operator, 
such as A, as a generator. Analogously to time evolu- 
tion, such a fiow has the form of a unitary operation, 
p{e) = exp(£(5^)(p(0)), where the real parameter e for 
the flow has units inverse to A. Therefore, taking the di- 
rectional derivative of pt along the specific curve of this 
fiow that passes through pi will produce ((39)) . For an 
explicit example that we will detail in ^VBl the state- 
space of a qubit can be parametrized as the continuous 
volume of points inside the unit Bloch sphere; the deriva- 
tive ([M)) produces the vector field illustrated in Figure [3] 
tangent to the flow corresponding to Rabi oscillations of 
the qubit. 

With this intuition in mind, we define the post- 
selection probability for measuring Pf given an initial 
state Pi{e) = exp{eSA){pi) that is changing along the 
fiow generated by A, 

Pf{e)=Tr,{PfP,{e)). (40) 



The logarithmic directional derivative of this post- 
selection probability then produces the factor 2ImA^ 
that appears in ([551). 

2W„ =5elnp/(£)|^^(j, (41) 

which is our main result. 

In words, the imaginary part of the weak value is half 
the logarithmic directional derivative of the post-selection 
probability along the natural unitary flow generated by 
A. It does not provide any information about the mea- 
surement of A as an observable, but rather indicates 
an instantaneous exponential rate of change in the post- 
selection probability due to disturbance of the initial state 
caused by A in its role as a generator for unitary trans- 
formations. Specifically, for small e we have the approx- 
imate relation, 

Pf{e)^Pf{0){l + {2lmA^)e). (42) 

For a pure initial state pi = \tlji){ipi\ and a projective 
post-selection Pf = \ipf){ipf\, the expression (pi]) simpli- 
fies, 

2ImA„ In 1(^/1 exp(-^el)|^,)|2|^^^J. (43) 

Hence, the corrections containing 2ImA^ that appear 
in (j38|) stem directly from how the specific von Neumann 
Hamiltonian ^ unitarily disturbs the initial system state 
infinitesimally prior to any additional disturbance in- 
duced by the measurement of the detector. Conceptu- 
ally, the coupling induces a natural unitary flow of the 
initial system state generated by A, which for infinitesi- 
mal g changes the joint probability p{p, f) for a specific 
p by the amount (2hiiAw){gp/h), where gp/h is the in- 
finitesimal parameter e in (|4ip that has units inverse to 
A. Averaging this correction to the joint probability with 
the detector observablcs i; or p produces the correction 
terms in ([55]). 



IV. TIME SYMMETRY 

As noted in [4^ [t^] , a quantum system that has been 
pre- and post-selected exhibits time symmetry. We can 
make the time symmetry more apparent in our opera- 
tional treatment by introducing the retrodictive state., 

Pf = P//Tr(P/), (44) 

associated with the post-selection (see, e.g., [2l,[zll) and 
rewriting our main results in the time-reversed retrodic- 
tive picture. 

After cancelling normalization factors, the detector re- 
sponse (j34p for a system retrodictively prepared in the 
final state pf that has been conditioned on the pre- 
selection measurement producing the initial system state 
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Pi has the form. 

Trs{£*{pf)p,) 
Tr,iV*iPf)p.) 



[XT 



9^^4MM^, (45a) 



2TT,i£*ipf)p,) 



i{p)t = 



TT,{£*{pf)p,) 



(45b) 



where the retrodictivc operations £*, X* , and V* are the 
adjoints of the predictive operations in ([27| and ([32]), 

t , 



(46a) 



X*{pf ) ^ {i,\{x,UlpfUT}l2\i,), (46b) 
V*{pf) = (V-Kp, [/;p/C/t}/2|^). (46e) 

Notably, the symmetric product with A that appears 
in the detector response (|45p involves the retrodictive 
state that has been evolved back to the initial time due 
to the non-selective measurement operation £* . Hence, 
in both pictures the measurement of A is being made 
with respect to the same initial time. 

After expanding the retrodictive operations perturba- 
tively as in 



ri=0 



(p")o(ad*A)"(p/), 



(47a) 



ri=0 



g \« ({p",a;})o 



(ad*A)"(/5/), (47b) 



(p"+i)o(ad*A)"(p/), (47c) 



where {&A* A){-) ~ — (adA)(-) = [•, A] is the right action 
of A, then the linear response of the detector (|38)) can be 
written in terms of the retrodictive forms of the real and 
imaginary parts of the complex weak value. 



Re-A^/; — 



2ImA„ 



Tr,({p/,A}p,) 
2Tr,(p/p,) ' 

Tvs{{~i[Pf,A])p,) 

9e In Tr^exp ( - e^A ) (p/ ) ) I _o , 



(48) 
(49) 



which should be compared with ([2|), dS]), and (|4T|) . 

We see that the imaginary part of the weak value can 
also be interpreted as half the logarithmic directional 
derivative of the pre-selection probability as the retro- 
dictive state changes in the opposite direction along the 
flow generated by A. 



V. EXAMPLES 

A. Bohmian Mechanics 

To make the preceding abstract discussion of the weak 
value more concrete, let us consider a special case that 



has been recently discussed by Leavens [58j, Wiseman 
[sgj . and Hiley [6l|, where the operator A = p being 
measured is the momentum operator of the system par- 
ticle. Since the wave-number operator k ~ —p/fi gener- 
ates a flow that is parametrized by the position a;, then 
we expect from the discussion surrounding (pij) that the 
imaginary part of a momentum weak value will give in- 
formation about how the post-selection probability will 
change along changes in position. 

If we restrict our initial system state to be a pure state 
Pi = 10) (01, and post-select the measurement of the mo- 
mentum on a particular position Pf ~ \x){x\, then the 
detector will have the linear response relations ([55)1 with 
the complex weak value given by. 



Pn 



~ihdx4>{x) 



(50) 



We can split this value naturally into its real and imagi- 
nary parts by considering the polar decomposition of the 
initial system state 0(a;) = r{x) exp(iS{x)), 



Pio = hdxS{x) — ihdx Inr(x). 



(51) 



The real part of the weak value Rep^, = hdxS{x) is 
the phase gradient, or Bohmian momentum for the ini- 
tial state, which we can now interpret operationally as the 
average momentum conditioned on the subsequent mea- 
surement of a particular x in the ideal limit of no mea- 
surement disturbance. This connection between the real 
part of a weak value and the Bohmian momentum that 
was pointed out in [sl, [H] has recently allowed Koc- 
sis et al. [sO] to experimentally reconstruct the averaged 
Bohmian trajectories in an optical two-slit interference 
experiment using such a von Neumann measurement. 

The imaginary part of the weak value, Imp^ = 
—hdx lnr(a;), on the other hand, is the logarithmic gradi- 
ent of the root of the probability density p{x) ~ |0(a;)P = 
r'^{x) for the particle at the point x. Written in the form. 



2Imp.u 



-hdx ln/9(a;), 



(52) 



it describes the instantanous exponential rate of posi- 
tional change of the probability density with respect to 
the particular post-selection point expected. This 

quantity, scaled by an inverse mass 1/to, was introduced 
under the name "osmotic velocity" in the context of a 
stochastic interpretation of quantum mechanics devel- 
oped by Nelson [s^l, where it produced a diffusion term 
in the stochastic equations of motion for a classical point 
particle with diffusion coefficient h/2m. Nelson's inter- 
pretation was carefully contrasted with a stochastic inter- 
pretation for the Bohmian pilot wave by Bohm and Hiley 
[8l| , and the connection of the osmotic velocity with a 
weak value was recently emphasized by Hiley [61| . 

Hence, the imaginary part of the momentum weak 
value does not provide information about a measurement 
of the momentum in the initial state. Instead, it indicates 
the logarithmic directional derivative of the probability 
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density for measuring x along the flow generated by p. 
The scaled derivative —hdx appears since p = —hk and 
k generates flow along the position x. 



Qubit Observable 



To make the full von Neumann measurement process 
more concrete, let us also consider a simple example 
where A operates on the two-dimensional Hilbert space 
of a qubit. (See also [H, M, M, We can 

in such a case simplify the perturbative expansions (|36p 
using the following identities, 



(53a) 

rko-kj , (53b) 

(53c) 
(53d) 



A 








pi 






= 2iejkio-i 







k 



where {(^fc}!^]^ are the usual Pauli operators, the com- 
ponents of the initial system state {rk}^^i arc real and 
satisfy the inequality < X^fe — ^jki is the com- 
pletely antisymmetric Levi-Civita pseudotensor, and Sjk 
is the Kroncckcr delta. Wc have deflned as to be diago- 
nal in the eigenbasis of A and have rescaled the spectrum 
of A for simplicity to zero out its maximally mixed mean 
Trs(Ai,72) = 0. As a result, {A)o = Ar^. 

It follows that for positive integer n the repeated ac- 
tions of A on the various qubit operators have the forms. 



(adA)"(i,) = 0, 
(adA)2"-i(6-i) = za-2(2^)2"-i, 

(adA)2»(^i) = ^i(2^)'", 
(adA)2"-i(^2) = -i^i(2A)2"-i, 

(adl)2"(^2) ^ ^2(2^1)2", 
(adA)"(^3) - 0, 



(54a) 
(54b) 
(54c) 
(54d) 
(54e) 
(54f) 



which collectively imply that, 



(adA)2"-i(pO = -(2A)2"-i (ri^2 - r^ai) 
(adA)2"(p,) = i(2A)2» (ri<Ti + r2&2) , 



(55a) 
(55b) 



has the exact form, 

- , c{g)ri - s{g)r2 



c{g) 

■<9) 



c{g)r2 + s{g)ri 



-(T2, 



^ (-1)" (2Ag^ 

n- 



oo 



E 



(-!)"+! {2Ag 



^ (2n-l)! 

Tl— 1 ^ ^ 



2n-l 



(56a) 

(56b) 

(56c) 



The correction term can be interpreted as a Rabi oscilla- 
tion of the qubit that has been perturbed by the coupling 
to the detector. Indeed, if the detector operator p were 
replaced with a constant p, then the interaction Hamilto- 
nian ^ would constitute an evolution term for the qubit 
that would induce Rabi oscillations around the axis of 
the Bloch sphere, which would be the natural flow in state 
space generated by the action of A. With the substitu- 
tion p —> p then (p")o — > p", so c{g) — > cos{2gAp/h) — 1 
and s{g) sm{2gAp/ti), which restores the unperturbed 
Rabi oscillations. 

Similarly, we find that the averaging operations for 
the detector position and momentum p2p have the exact 
forms, 



sx{g) 



Cx{g)ri - Sx{g)r2 



-(Tl 



(57a) 



Cxig)r2 + Sxig)ri 



^ i-ir ^2Ag^'" ({p2",4)o 



^ (2n)! V h 

^ (-1)"+^ ( 2Ag 



(57b) 



2n~l 



({p2»-l,4)o 



and, 



^ / \ - , Cp{g)ri -Sp{g)r2 



Cp{g)r2 + Sp{g)ri 



2ri 



cp{g) 
spig) 



^^^(2n-l)! V h J 



(57c) 
(58a) 

(58b) 
(58c) 



and hence that the nonselective measurement operation 



These operations differ from £ only in how the various 
moments of the initial detector distribution weight the se- 
ries for the Rabi oscillation. In particular, given the sub- 
stitutions p p and X ^ X, then ({p",x}/2)o — >■ p"a; 
and (p"+^)o p"+\ so Cx{g) x {cos{2gAp/h) - 1), 
Sxig) X sm{2gAp/h), Cp{g) p{cos{2gAp/h) - 1), 
and Sp{g) — > p sm{2gAp/h). Therefore, if the detec- 
tor remained uncorrelatcd with the system the averag- 
ing operations would reduce to X{pi) — > x£{pi) and 
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V{pi) — > p£{pi), which are the decoupled intial detec- 
tor means scaling the Rabi-oscillating qubit state. 

Since we have assumed that A does not have a compo- 
nent proportional to the identity, the symmetric product 
{A,pi}/2 = Ar3is/2 = (A)o(is/2) for a qubit wiU act 
effectively as an inner product that extracts the part of 
the initial state proportional to A. Therefore, the cor- 
rection to {A)o in that appears in (|34ap has the 
simple form, 



g 



TTs{Pf£{{A,p,}/2)) ^g{A)o 

Pif) ' 



(59) 



(60) 



where the conditioning factor, 

P{f)=2TrApf£{p^)), 
= 2Tr,(p/p,) 

+ {c{g)ri - s{g)r2)TT,{pf<Ti) 
+ (c(g)r2 -I- s{g)ri)TTs{pf(J2) 



is {2/TTs{Pf)) times the total probability of obtaining the 
post-selection. We have expressed p{f) more compactly 
in terms of the retrodictive state pi)) to show how the 
deviations from the initial state that are induced by A 
become effectively averaged by the post-selection state. 
In the absence of post-selection, the retrodictive state will 
be maximally mixed pf = 1^/2 andp(/) — > 1, recovering 
the unconditioned average {A)o. 

The correction to the detector mean position {x)q in 
(|34ap can be expressed in a similar way. 



Tr,(PjA-(pO) 1 /„, . 



Pif) 

+ {cx{9)ri - Sx{g)r2)Ti-s{pf&i) 
+ {cx{g)r2 + Sx(g)ri)Trs(p/o-2) 



(61) 



as can the correction to the detector mean momentum 
{p)o in 



Trs{PfV{p^)) 



1 



2{p)oTrs{pfPi 



+ {cp{g)ri - Sp(.g)r2)Trs(p/<Ti) 
+ {Cp{g)r2 + Sp(.g)ri)Trs(p/o-2) 



(62) 



Expanding ((Ml) using dSH), (HI]), and ^ to hnear 
order in 5, we find the linear response psp in terms of 
the real and imaginary parts of the qubit weak value. 



R,ej47M — 



2ImA„, = 



2Tr,(p/pO' 

Tv,{pf8A{p^)) 



^^^siPfPi) 

5a{,Pi) = A{-r2cri + 7'io-2). 



(63a) 

(63b) 
(63c) 



As expected, the real part contains information re- 
garding the measurement of A as an observable in the 




FIG. 3. (color online) The projection onto the plane rs = 
of the qubit Bloch sphere, showing the vector field 5a^[pi) = 
— r26"i + r\&2 for arbitrary initial states pi = (1 + ri<Ti + 
r2(T2 ~f r-iff-i)/2. The curves of the flow through this vector 
fleld are the Rabi oscillations around the axis that are 
generated by the unitary action of (T3. The quantity 2\raAw 
(|63p is the logarithmic rate of change of the post-selection 
probability (|41|) along this vector field. 



initial state, conditioned by the disturbance-free overlap 
between the predictive and retrodictive states. The imag- 
inary part, on the other hand, contains (^^(Pi), which is 
a tangent vector field on the Bloch sphere — illustrated in 
Figure [3] — that corresponds to an infinitesimal portion of 
the Rabi oscillation being generated by A. This tangent 
vector field contains only the components ri and r2 from 
bases orthogonal to A in the initial state pi, so 2ImA^ 
contains only the retrodictive averages of corrections to 
bases orthogonal to A, and thus contains no information 
about the measurement of A as an observable. As dis- 
cussed in (HH), 21mAyj is the logarithmic rate of change 
of the post-selection probability along the vector field 
^AiPi)- Scaling it by a small factor with units inverse 
to A will produce a probability correction to linear or- 
der. In the absence of post-selection, then pj ls/2, 
Re Aw ^ {A)o, and ImA^ 0. 



C. Gaussian Detector 

We can also apply our general results to the traditional 
case when the initial detector state in (fTO|) is a zero-mean 
Gaussian in position, 



(2^(T2)-i/4exp(-xV4a2), 



(64) 
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Then the measurement operators for position detection 
have the initial Gaussian form shifted by gA, 



1 



(271^2)1/4 



(65) 



while the conjugate measurement operators for momen- 
tum detection ((20|) have the initial Gaussian modified by 
a unitary factor containing A, 



The Wigner quasi-probability operator (|23ap corre- 
spondingly decouples into a product of Gaussian distri- 
butions, with only the position shifted by the system op- 
erator, 



(67) 



Marginalizing the Wigner operator over momentum 
and position separately produces the probability oper- 
ators (dH) and (PT|) . 



V27rcr2 



(68a) 
(68b) 



As anticipated, the probability operator for momentum 
no longer contains any information about the system op- 
erator A and is proportional to the identity, so measuring 
the momentum provides zero information about any sys- 
tem operator not proportional to the identity. 

In the presence of post-selection we can also exactly 
compute the disturbance operations (|27p . pOap . and 
(j30bp using the following identities for the Gaussian de- 
tector moments. 



(P'")0=(^) (2n-l)!!, (69a) 



(p'"-')o = 0, 
({p",x}/2)o = 0, 
(2n-l)!! 1 



(2n)! 



2"n!' 



(69b) 
(69c) 

(69d) 



which hold for positive integer n. We find the simple 
results, 

SiM = cxp [-i (A)' (adD^") (p,), (70a) 



7'(A) = |:^(adA)(f(pO). 



(70b) 
(70c) 



The quantity e = {g/2a)'^ with units inverse to A^ 
emerges as the natural decoherence parameter, which we 



can see more clearly by rewriting the non-selective mea- 
surement operation in (|70p as, 

/5,(e) = £{p,) = exp (eC[A]) (p,), (71a) 
C[A]{p{) = Ap,A^^]^{p,,A^A}. (71b) 

The operation C[A\{pi) is the Lindhlad operation [75l.[76l. 
HH that produces decoherence in continuous dynamical 
systems, with A playing the role of the Lindblad operator 
that decoheres the system. Since dtPi{e) = C[A]{pi{€)), 
the Gaussian measurement acts as an effective Lindblad 
evolution that decoheres the system state with increas- 
ing e via the action of A, but does not cause unitary 
disturbance along the natural flow of A [ssj . 

The exact expressions for the conditioned Gaussian de- 
tector means follow from ([M)) and ([70)1 . 



f{p)T 



Tr,(P^g({A,pJ)) 

2Tr,(Pff(p,)) ' 
g Tr,(Py(adA)(g(p,))) 



(72a) 
(72b) 



Surprisingly, the special properties of the Gaussian mo- 
ments (j69p allow (|72p to be written in a form proportional 
to the real and imaginary part of a complex weak-value 
involving the decohered system state ((7T|) to all orders in 
the coupling strength g, 



TrsiPfAp^ie)) 



Tr.(P/P,;(e)) 
f{x)T = gReAwie), 

/(p)T-|^(2W.(e)). 



(73a) 
(73b) 
(73c) 



Following the interpretations outlined in this paper we 
can therefore understand the position shift IleAyj{e) to 
all orders in g as the average of the observable A in the 
decohered initial system state pi{e) conditioned on the 
post-selection Pf. Similarly, we can understand the fac- 
tor 2ImA^(e) in the momentum shift to all orders in g as 
the logarithmic directional derivative of the probability of 
post-selecting Pf given the decohered initial system state 

Pi{e) along the unitary flow generated by A. 

If the measured operator is the qubit operator A = 
ASr^ as in (|53p . then we can further simplify the expres- 
sion ([70)1 using the identities ([5S)) to find, 

£{P^) =P^ + i(e-(^s/-)'/2 - 1) (ri,Ti + r2^2) , (74) 

= i (i + r3<T3 + e-(^^'/'^)'/2 (ri<Ti + r2<T2)) , 

which shows how the measurement decoheres the bases 
orthogonal to A in the initial state with an increase in the 
dimensionless flow parameter (Ag/cr)^ [g^. This deco- 
herence is illustrated in Figured The conditioned means 



13 




<5^,(p,(e)) 




FIG. 4. (color online) Two projections of the Bloch sphere 
showing the pure decoherence of the specific state pi[e) — 
exp(e>C[<T3])(pi) = (i + exp(^2e)x/3o-2/2 + d-3/2)/2 due to 
the Gaussian detector (|74p . (left) The projection onto the 
plane ri = showing the progressive collapse of pi{e) onto 
the rs axis with increasing e. (right) The projection onto the 
plane rs = showing the vector field Scr^ipii^)) during the 
progressive collapse. Notably the quantity 2Imylu,(e) (|73p is 
the rate of change of the post-selection probability (141(1 along 
this vector field for all e, but not along the purely decohering 
trajectory that Pi{e) actually follows. 



(17^ of a Gaussian qubit detector consequently have the 
exact form, 



f{x)T 



'p(/)' 



(75a) 
(75b) 



(riTr,(p/o-2) - r2Tr^(p/<Ti)), 
p(/) = l + r3Tr,(p/^3) (75c) 



+ e' 



(riTr,(p/0-i) + r2Tr,(p/^2)) 



to all orders in the coupling strength g. When expanded 
to linear order in g, (|75|) reduces to (|38p with the real 
and imaginary parts of the qubit weak value (j63p . as 
expected. 

For contrast, as g becomes large the unconditioned 
measurement of A becomes essentially projective and the 
operation £ almost completely decoheres the initial state 
([7^ into the basis of A as the pointer basis, 



^(pO ~ 2 (is +^3^3) 



(76) 



Hence, in this strong measurement regime, the condi- 
tioned means (|75p approximate. 



(77a) 
(77b) 



^ {A)o 
^<^^---^l + .3Tr.(p,^3)' 



0. 



The position shift contains the average of A in the deco- 
hered initial system state £{pi), conditioned by the post- 
selection. Moreover, since the decohered initial system 



state £{pi) is essentially diagonal in the basis of A, it 
will no longer Rabi oscillate, so the directional derivative 
along the flow generated by A will be essentially zero. 
Hence, the probability correction factor represented by 
2ImA„(e) vanishes. 



VI. CONCLUSION 



We have given an exact treatment of a conditioned von 
Neumann measurement for an arbitrary initial state and 
an arbitrary post-selection using the language of quan- 
tum operations. The full form of the conditioned de- 
tector response (p4|) naturally indicates how the mea- 
surement disturbance and conditioning from the post- 
selection modify the unconditioned detector response. 
The corresponding linear response of the detector ([38)) 
can be parametrized by the generalized complex weak 
value ([35]) , but the origins of the real and imaginary parts 
differ. 

The real part of the weak value ([2]) stems directly from 
the measurement of A made on the initial system state, 
and can be interpreted as the idealized zero-disturbance 
conditioned average of the operator A acting in its role 
as an observable. The imaginary part of the weak value 
Q ; on the other hand, contains no information about the 
measurement of A as an observable, but instead arises 
from the disturbance due to the von Neumann coupling. 
We interpret it as the logarithmic directional derivative 
(|4T|) of the post-selection probability along the unitary 
flow in state space generated by the operator A in its 
role as the element of a Lie algebra. The complex weak 
value therefore captures both halves of the dual role of 
the operator A in the quantum formalism. 

To illustrate this interpretation, we considered the 
weak value of momentum (j51|) post-sclcctcd on a par- 
ticular position. Its real part is the Bohmian momen- 
tum representing the average momentum conditioned on 
a position detection, while its imaginary part (|52p is pro- 
portional to the "osmotic velocity" that describes the 
logarithmic derivative of the probability density for mea- 
suring the particular position directed along the flow gen- 
erated by the momentum. 

Finally, we applied our exact solution to the useful 
special cases of a qubit operator and a Gaussian detec- 
tor. Wc showed how the natural qubit Rabi oscillations 
that would be generated by the von Neumann interaction 
Hamiltonian (j4|) become disrupted by the measurement. 
We also showed that the Gaussian detector purely de- 
coheres the initial system state into the basis of A in 
the Lindblad sense, which allows the exact interpreta- 
tion of the position and momentum shifts of the detector 
in terms of a complex weak value (|73p that involves the 
decohered system state. 
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